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Abstract 

Recent theoretical results on heavy flavor production and 
decay in the framework of perturbative QCD are reviewed. 
This includes calculations for top production at hadron collid- 
ers, inclusive charmonium production and the comparison be- 
tween the singlet and octet mechanisms. Predictions for heavy 
flavor production in e^e~ annihilation will be discussed in some 
detail, covering both the threshold and the high energy region. 
The first results in NLO for heavy flavor decays will also be 
reviewed. 



The complete paper, including fig ures, is also available via anonymous ftp at 



ftp: / / ttpux2.physik.uni-karlsruhe.de/ 



or via www at 



http: / /www-ttp. physik.uni-karlsruhe.de / cgi-bin/preprints / 



■^Supported by BMBF contract 057KA92P. 
Invited talk presented at the "Cracow International Symposium on Radiative Cor- 
rections", Cracow, Poland, 1996. 



1 Introduction 



Heavy flavor production and decay have developed into benchmark 
reactions for perturbative QCD. The large energy scale inherent in 
most of these reactions allows for a separation between hard and soft 
momentum transfers. The former can be treated perturbatively, the 
nonperturbative matrix elements which encode the remaining infor- 
mation can either be determined experimentally, or integrated out by 
considering sufficiently inclusive information such that perturbation 
theory alone is adequate. 

Significant progress has been achieved recently in a number of top- 
ics. The predictions for top production at hadron colliders have been 
scrutinized by several authors. In particular the role of soft gluon 
resummation has been emphasised and the ag dependence explored 
(Section 2). Inclusive charmonium production at hadron and e+e" 
colliders has been studied theoretically and experimentally. A fairly 
complex picture seems to emerge, with different mechanisms playing 
a role in various reactions (Section 3). The inclusive cross section for 
heavy fiavor production in e^e~ annihilation has been studied in a va- 
riety of papers. Far above threshold an expansion in / s is adequate 
and has been successfully applied to Z decays to bottom quarks, or to 
charm production just below the bb threshold. For a prediction above, 
but relatively close to threshold a different strategy has been employed, 
which is based on a combination of analytical and numerical methods. 
For an adequate treatment of top quark production in the threshold 
region its large decay rate and the interplay between gluon radiation 
from the production and the decay process must be taken into ac- 
count. These topics will be reviewed in section 4. The leading QCD 
correction to weak decays of heavy fiavors have been evaluated quite 
some time ago. Results are available for the rate, the spectrum and 
for angular distributions. To match the level of precision claimed by 
the proponents of the Heavy Quark Effective Theory, next to leading 
order predictions are required from perturbation theory. First steps 
into this direction have been made and will be reviewed in section 5. 
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2 Top production in hadronic collisions 



The theoretical framework and the (semi-) analytical results for the 
top production cross section in NLO have been developed nearly a 
decade ago ^ 0. The predictions for y/s =1.8 GeV and rrit = 180 
GeV from various authors are listed in Table |l]. For = 175 GeV the 
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Table 1: History of predictions for the production cross section for 
= 1.8 TeV and = 180 GeV. 

cross section increases by about 0.7 pb. The uncertainty in the factor- 
ization and renormalisation scale leads to an uncertainty of roughly 
10%. Recently the issue of soft gluon resummation has been raised. 
The original arguments ^] leading to a large positive shift of roughly 
10% have been refuted in 0. No consensus has yet been reached on 
the magnitude of these effects. Increasing from the nominal value 
of around 0.11, which has been frequently used in these calculations, 
to 0.120 leads to an increase by about 5%. Within the combined un- 
certainties theory and experiment are in very good agreement (Fig. |l]). 

All these calculations are based on a perturbative treatment of 



3 




4 



the threshold region. In principle one should, however, incorporate 
the leading terms of order TTas/j3. The resulting modifications are 
small for ti in a color octet which is the dominant configuration at the 
TEVATRON (see section 2.2.2 in H). 



3 Inclusive Charmonium Production 

High energy hadron-hadron and e~p colliders are charmonium facto- 
ries. A variety of production mechanisms have been discussed in the 
literature. Contributing with different relative strengths in the various 
reactions they can be disentangled only through a systematic study of 
different processes. In particular the question of color singlet versus 
octet production has stimulated a number of detailed investigations. 

Inelastic J /ip production in photon-photon reactions provides a rel- 
atively clean testing ground. The dominant subprocess at the parton 
level 

-i + g^J/i^ + g (1) 

can produce directly a (cc) color singlet state. Incorporating also the 
one loop perturbative corrections 0, satisfactory agreement between 
theory and experiment is observed for the J/ip energy distributions 
and the total production cross section as well (Fig. H). 

This success of the color singlet model (CSM) (where quarko- 
nium (color singlet!) states are required to be produced through a 
purely perturbative mechanism) is in marked contrast with its failure 
in purely hadronic collisions. The dominant subprocesses in the CSM 
are based on the conversion of a virtual gluon into J/ if) or xj plus 
two or one gluon respectively. The combination of additional pow- 
ers of as with the small phase space gives rise to sizable suppression 
factors. This perturbative treatment of soft gluon radiation may be 
inadequate and an alternative approach has been advocated in . 



The cross section for charmonium production is decomposed into a 
sum of terms consisting of the cross section for (cc) states in a specific 
angular momentum and color state times the nonperturbative matrix 
element of an operator characterizing the conversion probability into 
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Figure 2: Comparison between theoretical prediction for the energy 
dependence of the inelastic cross section {z < 0.9) for J/ip photo- 
production (J. Steegborn, private communication, based on P) and 
recent data from the HI Collaboration. 

J/^: 

a{pp ^ J/^ + x) = Y^ a{pp ^ cc{n)) x (C>;[/^). (2) 

n 

These matrix elements are effectively free parameters to be deter- 
mined in different experiments. This approach is closely related in 
its spirit to the color evaporation model formulated a long time ago; 
it provides, however, a more firm theoretical formulation. Adjusting 
the parameters appropriately, a satisfactory description of the data is 
obtained. 

The clean initial state configuration typical for e+e~ annihilation 
is ideal to investigate the relative importance of different production 
mechanisms. Two distinctly different situations have been considered: 
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high energy reactions hke Z decays with large event rates available 
at LEP and alternatively the 10 GeV region that can be explored 
at present at CESR or in the near future at the B-meson factories. 
Three mechanisms have been identified at which contribute in the high 



energy region with comparable rates. The reaction |TT 



Z J/iPcc + X (3) 

requires the production of two cc pairs with a rate proportional to 
ag|-R(0)p. The second mechanism ||12| is the splitting of a virtual 
gluon in a color octet cc: 

Z qq{cc)s (4) 

with the subsequent nonperturbative conversion of (cc)8 into J/ip. The 
rate for this mechanism is proportional to where the second 

factor characterizes the nonperturbative matrix element. The third, 
color singlet, contribution |jT3| 



Z^qqJ/ipgg (5) 

is strongly suppressed by the factor ag|i?(0)p and, furthermore, by 
the small phase space. The branching ratios of the three reactions 
are given by 0.8 ■ 10-^ 1.9 ■ 10~^, 0.2 ■ 10"^, respectively. The total 
inclusive rate is reasonably consistent with the observations by the 



OPAL collaboration flj of (1.9 ± 0.7 ± 0.5 ± 0.5) ■ 10"^ However, 
a statement about the dominance of any of these processes seems 
premature. The analysis of J/ip energy and momentum distributions, 
however, could help to settle this issue. 

Also B meson factories and CESR give rise to a large sample of 
events with prompt J/ ip production. Two mechanisms have been pro- 
posed which might well describe complementary kinematical regions. 
The leading process in the CSM 

e^e~^J/ip + gg (6) 

is proportional to ag|i?(0)p. It leads to a three body final state and 
hence to a continuous energy distribution (Fig. 0). Predictions for the 
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Figure 3: Energy distribution for inclusive J /if) production in e^e 
annihilation at 10.5 GeV. Solid curve: q;s(M| ), dashed curve as^Ml). 



rate, the angular and the momentum distribution and the polarization 
can be found in [|T^. The alternative approach is based on "color 
octet production", e+e~ — ^ {cc)% + g. The rate is of order as and 
multiplied by a nonperturbative matrix element. The J/ip energy is 
essentially fixed at Emax = (-^ + "^^)/(2A/i). The angular distribution 
is proportional to (1 + cos^0). These features are identical to the 



predictions of the "color evaporation model" IjT^ formulated a long 
time ago. An excess of J /ip at this special kinematical point with the 
predicted angular distribution would be a strong indication for this 
"octet mechanism" . The angular distribution of the J/ in the GSM 
is of the form 1 + a{y) cos^ 9 where a{y) depends on y = E^/EBeam 
and approaches roughly —0.8 at the endpoint (Fig. ^). This difference 
will be crucial in disentangling the two mechanisms. 
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Figure 4: Coefficient a{y) characterizing the angular distribution of 
J/tjj^s (L: longitudinally polarized J/'^'s only). 



4 Heavy flavor production in e+e anni- 
hilation 

4.1 Z ^bb 

Experimental studies of various partial and of the total Z decay rate 
have been performed recently with a new level of sophistication. The 
relative error in Tb has been lowered to about 0.5 ■ 10~^ corresponding 
to STf) ^2.5 MeV, the uncertainty in the total decay rate which is also 
influenced by F;, amounts to about 3 MeV. In comparison with F^ or 
Fy two important differences have to be taken into account for F;,. 
The first, relatively straightforward aspect is related to the bottom 
mass. In Born approximation the correction from the phase space 
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suppression of the axial part of the rate is predicted to be — 6m^/M| 



corresponding to —4 MeV. In |jT8[ it has been demonstrated that this 
number is drastically modified by QCD corrections. The bulk of these, 
the large logarithms, can be absorbed by reexpressing the result in 
terms of the running mass thus reducing the correction to —1.6 MeV. 
(For a detailed discussion and further references see |1^.) The second 
contribution to the Z —>■ bb decay has its origin in the double triangle 
diagrams with two gluon intermediate states. It is present for the 
axial rate only. The contribution of order was calculated quite 
some time ago for arbitrary mf/M'^. Formally it is proportional to 
lnm^/M| and thus seems to diverge in the limit of large In m^/M|. 
However, additional logarithms of m^//!^ are induced by the running 
of as which have to be controlled at the same time. The structure 
of leading logs was analysed in pD|, the constant terms of in El 



The combined effect of order and from these "singlet terms" 
amounts STf, = — 1.8 MeV. It is clear that the sum of mass and singlet 
terms must be taken into consideration in any precision analysis. 



4.2 Intermediate energies 

The Z decay rate is well described in the massless approximation plus 
terms of order ml/M'^. However, for a prediction at lower energies, 
an increasing series of terms in the m^/s expansion is needed. The 
comparison between the complete calculation and a limited number 
of terms in the m^/s expansion indicates that the first three terms 
are sufficient to describe the cross section from high energies down to 
s ~ 8m^. With this motivation in mind the quartic terms of order 

have been calculated in In this way an adequate prediction 
between roughly 14 GeV and Mz is available for bb production, and 
similarly for cc production from roughly 5 to 6 GeV up to the bottom 
quark threshold ||23| (Fig. |^). 

In view of the large statistics available at CESR and at a future 
5-meson factory a detailed theoretical study has been performed in 
[ P^ which demonstrates the potential for this potentially most precise 
and clean determination of a,. 
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Figure 5: The ratio R{s) below and above the b quark production 
threshold at 10.5 GeV for a,(Mz) = 0.120, 0.125 and 0.130. The 
contributions from light quarks {u, d, s, c) and the bottom quark are 
displayed separately. 



4.3 The NLO calculation for arbitrary and s 

A few GeV above charm, bottom, or top threshold measurements can 
in principle be performed at a r-charm factory, at a i?— meson factory 
and a future linear collider. With a relative momentum of the quarks 
exceeding for instance 3 GeV perturbative QCD should be applicable 
also in this region. It is, therefore, desirable to push the theoretical 
prediction as close as possible towards the threshold. The two-loop 
calculation has been performed more than 40 years ago . The imag- 
inary part of those three-loop diagrams which originate from massless 
fermion loop insertions in the gluon propagator ("double bubble dia- 
grams") were calculated analytically in |2^. Real and imaginary parts 



of the purely gluonic correction (and of the double bubble diagrams) 
were calculated in a semianalytical approach that will be sketched 
in this subsection. 
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The polarization function can be written in the form 

n = n(o) + ^n(i) 

vr 

+ (^) (c|nA + CfCaHna + CpTnilii + CfTHf) (7) 

where ni denotes the number of hght quark species. Each one of the 
Ilj is analytical in the complex plane with a cut from Am? to +00. 
For small they can be expanded in a Taylor series 

The renormalization condition n(g^ = 0, m?) = has already been 
implemented. The evaluation of the Taylor coefficients amounts to the 
calculation of three loop tadpole integrals with an increasing number 
of mass insertions - up to 16 for Cg which is the present limit for the 
evaluation with the help of algebraic programs. 

In the large region a similar expansion can be performed. For 
this case the expansion has been performed up to terms of order 
{w? /q^)^ and {m? /q^Y- Additional information can be obtained about 
the behavior close to threshold. Leading and subleading terms can be 
deduced from the influence of the Coulomb potential in the nonrel- 
ativistic region, combined with the knowledge about the logarithmic 
corrections of the perturbative QCD potential. To extend the range 
of convergence from q^ < Am? to the full analyticity domain an ap- 
propriate variable transformation has to be performed. The data from 
the three kinematical regions are flnally integrated in a Fade approx- 
imation which leads to stable results for n(g^) and R{s) at the same 
time. The result for the three dominant pieces are shown in Fig. ^ 
where it is compared to the leading terms close to the threshold and 
to the high energy approximation. 

4.4 Toponium and top quarks in the threshold 
region 

Top quarks were treated as stable particles in the previous section. 
Although adequate away from the threshold, this approximation is in- 
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Figure 6: Complete results plotted against v = yl — im^/s. The 
high energy approximation includes the rn^ / term. 
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adequate in the "would-be" toponium region. For a mass of the top 
quark around 175 GeV a decay rate Tt ~ 1.5 GeV is predicted, cor- 
responding to a toponium width of 3 GeV. The resonances are thus 



completely dissolved p8|, and the individual peaks are merged 
into a step function like threshold cross section. Quarkonium physics 
ceases to exist. The large decay rate introduces, however, a cutoff 
which eliminates all nonperturbative aspects of the interquark poten- 
tial. Large momentum tails beyond 



-cut - \/2mtTt ^ 24GeV (9) 
or, alternatively, distances above 

r ^ O.Olfm (10) 



are irrelevant for the description of the tt system ||3y, 0, |33|. The 
impact of the large rate is clearly visible in Fig. |^. The predictions 
for three different top masses rrit = 150 GeV, 180 GeV, and 200 GeV 
corresponding to = 0.81 GeV, 1.57 GeV, and 2.24 GeV demonstrate 
the strong influence of F^ on the shape of the cross section. The shape 
is furthermore significantly modified by initial state radiation and the 
spread in the beam energy. 

Additional information is encoded in the momentum distribution 
of top quarks, the "Fermi motion" which can be traced through the 
decay products W + b. This distribution is essentially equivalent to the 
square of the wave function in momentum space and can, for unstable 
particles, be evaluated ^ with Green's function techniques 

(Fig. H). Various experimental studies have demonstrated the poten- 
tial of a linear collider to determine rrit to a precision of perhaps even 
200 MeV by measuring the cross section and the momentum distribu- 
tion simultaneously. 

Highly polarized top quarks are required for a variety of precision 
studies of top decays. In the threshold region this is easily achieved. 
In fact, even with unpolarized beams top quarks are longitudinally 
polarized (with a polarization around —0.4) as a consequence of the 
nonvanishing axial part of the neutral current. Longitudinally polar- 
ized beams lead to a fully polarized sample of top quarks. 
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Figure 7: Total cross section as function of E = ^/s — 2mt for three 
values of the top quark mass. 



Another step in complication is achieved by considering the inter- 
ference between the dominant S and the suppressed P wave contribu- 
tions. The relative size of these effects is of order /5 ~ 0.1. It leads 
to a forward-backward asymmetry and furthermore to an angu- 
lar dependent quark polarization perpendicular to the beam direction 
p5[] . A detailed discussion of these effects, in particular of the role of 
the normal polarization and of rescattering corrections, can be found 
in . The small polarization of top quarks normal to the production 
plane is a particularly sensitive measure of the interquark potential. 

Additional complications are introduced through the rescattering 
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Figure 8: Momentum distribution of top quarks for three different cms 
energies. 



[ P^ , P7[ between b quark jets and the spectator, and by relativistic 
corrections [^] of order a^. These effects will be important for the 



quantitative comparison between theory and experiment and the ex- 
traction of a precise value for rrit, Ff and from threshold studies. 



5 Towards NLO in Heavy Flavor Decays 

Semileptonic weak decays of bottom mesons and top quarks are par- 
ticularly clean probes of the fundamental properties of quarks, their 
masses and mixing angles. Decay rates are, however, influenced by 
QCD effects, a large part of which can be calculated in PQCD. Leading 
order corrections to practically all quantities of interest are available: 
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for the decay rate of charmed and bottom quark from and for top 
quarks from [^]. Lepton decay spectra have been calculated in |40|| , 
the energy distribution of hadrons in . Leptons from the decays of 
polarized quarks exhibit a nontrivial angular distribution ^ and 
even lepton mass effects have been incorporated in these calculations 
|45| . A compact summary of most of these QCD corrections can 
be found in |^ . 

Different techniques to determine the degree of b or top polarization 
have been investigated in The analysis of moments of the 

lepton momentum distribution, or the ratio of charged vs. neutral 
lepton moments appear to be particularly promising. 

The corrections are often sizable, in particular those to the decay 
rate. In order to fix the scale in the running coupling constant and to 
gain confidence in the numerical result, a calculation of NLO correc- 
tions to the rate, if not the spectrum, is necessary. Purely fermionic 
loops have been considered in [|^, |49(| . 
result is particularly simple 
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(12) 



If we adopt the BLM prescription the large coefficient leads to a large 
shift in the effective scale for a^: fisLM = 0.12mt. Similarly large 
correction factors have been observed ||5^ for the decay of b into lu 
plus a charmed or u quark. 

It should be emphasized that the magnitude of NLO corrections 
~ (as(mj)/7r)^ ~ (0.07)^ is well comparable with correction terms ob- 
tained in Heavy Quark Effective Theory — typically of order (A / rrih)'^ ^ 
(0.05)^. Transitions at zero recoil i.e. for the final state with pc = ^Pb, 
are particularly clean from the theoretical point of view. No uncalcu- 
lable form factor is present, allowing to determine Vcb with remarkable 
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precision. The first calculation of the full NLO QCD corrections has 
therefore been performed at zero recoil [51]. Two important simplifi- 
cations are present in this case: 



• no real radiation has to be considered, 

• only relatively simple two loop integrals arise which can be cal- 
culated in a series expansion. 

The resulting NLO corrections are smaller than the leading ones 
by about a factor 4, reducing thus the theoretical error by a significant 
factor. Evidently these results can be considered a first important step 
towards a complete NLO calculation of the heavy quark decay rate. 
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